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ABSTRACT 

We report the oxygen abundances of the H II regions of a sample of low surface brightness 
(LSB) galaxies. We provide analytic functions describing the McGaugh (1991) calibration 
of the i?23 method. We use this and the equivalent width (EW) method to determine oxygen 
abundances, and also make direct estimates in a few cases where the temperature sensitive 
[O III] A4363 line is available. We find LSB galaxies to be metal poor, consistent with the L-Z 
relation of other galaxies. The large gas mass fractions of these objects provide an interesting 
test of chemical evolution models. We find no obvious deviation from the closed-box model 
of galactic chemical evolution. Based on our abundance and gas mass fraction measurements, 
we infer that LSB galaxies are not fundamentally different than other galaxy types but are 
perhaps at an early stage of evolution. 
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1 INTRODUCTION 

Low surface brightness (LSB) galaxies are normal sized disk 
galaxies with peak surface brightnesses of Ho ^ 23B mag 
arcsec^^. The sky background makes it very difficult to study 
the weak signal from their stellar continuum. Emission lines from 
bright H II regions embedded within these galaxies provide more 
practical probes. Prominent among these are the strong lines of 
oxygen from which the abundance of that element can be deduced. 
As a primary product of Type II supernovae and the most abun- 
dant element after hydrogen and helium, oxygen provides a useful 
measure of the degree of chemical evolution of a galaxy. 

LSB galaxies tend to have gas fractions that are systematically 
larger than well- studied brig ht galaxies iMcGaugh & de Blokl 
ll997HSchombert. McGaugh. & Eder 200ll) . Since gas fraction is 
the primary variable in theories of chemical evolution, LSB galax- 
ies afford the opportunity to test models like closed-box evolu- 
tion in an important regime. This also provides a check on the 
luminosity-metallicity (L-Z) relation determined for other galaxy 
types. 

In this paper we study the optical emission line spectra from 
H II regions in a number of LSB galaxies. In Section 2 we discuss 
the observations and data reduction. In Section 3 we compare our 
spectra to diagnostic diagrams. Oxygen abundances are calculated 
in Section 4 and discussed in Section 5. The luminosity-metallicity 
and mass-metallicity relations are discussed in Section 6, followed 
by the gas mass fraction in Section 7. Conclusions are presented in 
Section 8. 



2 OBSERVATIONS AND DATA REDUCTION 

The William Herschel Telescope at La Palma was used to obtain 
optical spectra of H II regions in six low surface brightness (LSB) 
galaxies during the nights of 1995 January 4-6. A 2.4 arcsec slit 
with 7 A resolution was used to take blue spectra spanning 3650- 
5250 A and red spectra spanning 5700-7300 A. Five standard stars 
and CuNe and CuAr calibration lamps were also observed for flux 
and wavelength calibration, respectively. 

The LSB galaxies observed were selected for containing rea- 
sonably bright, low metallicity Hll regions (as determined by 
lMcGaughlll994) . The slit was placed along each galaxy so its po- 
sition would maximize the number of H II regions observed during 
each exposure. There were two slit positions for UGC 12695. Indi- 
vidual exposures were 900 seconds, and multiple exposures were 
taken to increase the signal-to-noise and reject cosmic rays. Cumu- 
lative exposure times ranged from 7200 s to 20,700 s. Table 1 lists 
the galaxies and number of H II regions observed. Images of the 
six galaxies observed are shown in Figs. 1-6. North is up and East 
is to the left in the images. H II regions with oxygen abundance 
measurements are labeled in each image. Further discussion of the 
locations of the Hll regions can be found in McGaugh (1992), 
and more images are published in McGaugh, Schombert & Bothun 
(1995). 

Standard data reduction techniques were carried out using 
IRAF ^. First, the average bias frame was subtracted. Next, sky 
flats were used to ensure uniform illumination along the slit and 
tungsten flats were used to flatten the observations in the spectral 



^ The Image Reduction and Analysis Facility (IRAF) is distributed by the 
* Email: kuzio@astro.umd.edu (RKD); ssm@astro.umd.edu (SSM); Er- Association of Universities for Research in Astronomy, Inc., under contract 

win.deBlok@astro.cf.ac.uk (EdB) to the National Science Foundation. 
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Table 1. Summary of Observations. 



Galaxy 


No. of H II Regions 


Total No. Exposures 


F563-1 


3 


23 


F571-5 


1 


9 


UGC 1230 


5 


20 


UGC 5005 


2 


14 


UGC 9024 


2 


19 


UGC 12695 (1) 


2 


8 


UGC 12695 (2) 


3 


9 



Two slit positions were used for UGC 12695. 
Individual exposures were 900 seconds. 




Figure 1. i?-band image of F563-1. North is up and East to the left. The 
observed H II regions are labeled along the slit. 




Figure 2. V-band image of F571-5. The observed H II region is labeled. 
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Figure 3. Ha image of UGC 1230. Observed H II regions are labeled along 
the slit. 
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Figure 4. J?-band image of UGC 5005. The observed H II regions are la- 
beled along the slit. 



direction. Wavelength calibration of the observations was achieved 
by identifying spectral lines in the CuNe and CuAr calibration 
spectra. The Hll region spectra were flux calibrated using the 
standard star observations. Lastly, the individual H II region spec- 
tra were extracted from the galaxy frames and saved as one- 
dimensional spectra. The spectra from multiple observations were 
co-added. Fig. 7 shows an example of an H II region spectrum with 
important nebular lines identified. It should be noted that this is one 
of the better spectra taken and is not typical of the entire sample. 

There is no overlap between the blue and red ends of the spec- 
tra. Based on previous experience with the instrument, we antici- 
pated no difficulty in reconciling the calibrations of the red and 
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Figure 5. Ha image of UGC 9024. Observed H II regions are labeled along 
the slit. 



Figure 7. Spectrum of HII region 2 in UGC 1 2695 slit position 1 . Important 
nebular lines are labeled. 
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ic ■ . ■ ■ ■ 

Figure 6. Ha image of UGC 12695. Observed Hll regions for both slit 
positions are labeled. Slit position (1) runs East-West; slit position (2) runs 
roughly North-South. 



blue spectra with standard reduction procedures. While this is of- 
ten the case, we found that the flux scales are not always consistent 
in that they have a mismatch between continuum levels or give 
unphysically low Hq/H/3 ratios (< 2). This difference between 
arms is inconsistent among different galaxies and sometimes even 
between Hll regions of the same galaxy. We therefore have no 
confidence in the absolute flux calibration of the data, nor can 
we confidently connect the flux scales of the red and blue spec- 
tra. However, the line ratios within each spectrum (e.g., [O IIl]/H/3 
and [N Il]/Ha) can be measured even if those between red and blue 
spectra (e.g., Hq/H/9 and [N Il]/[0 II]) cannot. We therefore restrict 
our analysis to line ratios and equivalent widths from within each 



spectrum. The important abundance sensitive lines are all in the 
blue spectra, so this is not a great impediment. The most serious 
problem which arises is that the reddening cannot be estimated 
from the Ha/H/3 ratio, though in many cases single spectra con- 
taining both these lines already exist (McGaugh 1994). We list in 
Table 2 the observed blue line fluxes relative to H/3 and in Table 3, 
the observed red line fluxes relative to Ha. 



3 DIAGNOSTIC DIAGRAMS 

Considering the difficulties with reconciling the flux scale between 
red and blue spectra, it is worth making some checks to confirm 
that line ratios within each of the red and blue spectra do in fact 
make sense. This can be done by comparison to independent data 
(McGaugh 1994) and with diagnostic diagrams a la Veilleux & 
Osterbrock (1987) (Fig. 8). The line in Fig. 8 is the relation for H II 
regions from Baldwin, Phillips & Terlevich (1981). In this plot, 
AGN fall in the upper right and H II regions along the line. Fig. 8 
confirms that our spectra have line ratios consistent with those of 
H II regions. 

Having confirmed that we have H II region spectra, we now 
compare the current data to the data of McGaugh (1994). In Fig. 9 
we plot the ratio of the [O II] and [O III] lines of the current sample 
to the McGaugh (1994) sample against the flux of the H/3 line from 
McGaugh (1994). If the same regions have been observed then the 
ratios should be the same. It is clear from the figure that, with the 
exception of two discrepant H II regions in the [O II] line, the same 



© year RAS, MNRAS OOO.ITIfni 



Table 2. Observed "Blue" Emission Line Intensities Relative to Hp. 



Galaxy 


E563-1 


F563-1 


F563-1 


F571-5 


U1230 


U1230 


U1230 


U1230 


U1230 


Hll 


1 


2 


3 


1 


1 


2 


3 


5 


6 


[Oll]A3727 


4.15±0.29 


2.36±0.17 


5.63±0.40 


2.41±0.17 


1.71±0.12 


2.58±0.18 


1.76±0.12 


1.61±0.11 


3.78±0.27 


[Neiii]A3869 




0.45±0.03 


0.78±0.06 


0.31±0.02 




0.23±0.02 


0.28±0.02 






Hi;A3889 




0.19±0.01 


0.39±0.03 








O.lOiO.Ol 






HeA3970 




0.20±0.01 


1.09±0.08 








O.lOiO.Ol 






[Sii]A4076 




















H5A4102 




0.23±0.02 


1.43±0.10 


0.44±0.03 


0.29±0.02 


0.25±0.02 


0.21±0.02 


0.64±0.05 


0.82±0.06 


H7A4340 


0.51±0.04 


0.48±0.03 


1.55±0.11 


0.12±0.01 


0.34±0.02 


0.43±0.03 


0.33±0.02 






[Oiii]A4363 


0.16±0.01 


















[Oiii]A4959 


0.69±0.05 


1.81±0.13 


1.00±0.07 


1.11±0.08 


0.85±0.06 


0.36±0.03 


0.87±0.06 


0.83±0.06 


0.25±0.02 


[Oiii]A5007 


1.87±0.13 


5.22±0.37 


2.44±0.17 


3.48±0.25 


2.21±0.16 


1.54±0.11 


2.70±0.19 


2.77±0.20 


3.59±0.25 


F(H/3) 


0.71±0.04 


3.11±0.16 


0.26±0.01 


1.43±0.07 


0.37±0.02 


1.17±0.06 


2.23±0.11 


0.27±0.01 


0.16±0.01 


EW(H^) 


36±2 


21±1 


13±1 


16±1 


59±3 


41±2 


49±2 


13±1 


6±1 


EW([Oll] A3727) 


140±7 


69±3 


87±4 


36±2 


67±3 


131±6 


lOliS 


31±2 


30±2 


EW([OlIl] A4959) 


27±1 


38±1 


13±1 


15±1 


39±2 


14±1 


35±2 


9±1 


1±1 


EW([Olll] A5007) 


76±4 


111±6 


31±2 


54±3 


102±5 


56±3 


105±5 


33±2 


20±1 


H/3 flux in 10-1= , 


erg cm^^ sec" 


1; EW in A. 
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Table 2 - continued 



@ 



Galaxy 


U5005 


U5005 


U9024 


U9024 


U12695(l) 


U12695(l) 


U12695(2) 


U12695(2) 


U12695(2) 


Hii 


2 


3 


1 


3 


2 


3 


2 


3 


4 


[Oii]A3727 


3.38±0.24 


2.74±0.19 


3.50±0.25 


2.00±0.14 


2.28±0.16 


3.96±0.28 


1.78±0.13 


1.49±0.11 


1.15±0.08 


[NeIIl]A3869 


0.31±0.02 


0.33±0.02 




0.43±0.03 


0.35±0.03 


0.83±0.06 


0.39±0.03 


0.24±0.02 


0.28±0.02 


HCA3889 


0.21±0.02 


0.17±0.01 






0.14±0.01 


0.24±0.02 


0.17±0.01 


O.lSiO.Ol 


0.16±0.01 


HeA3970 


0.22±0.02 


0.19±0.01 


0.27±0.02 


0.22±0.02 


0.25±0.02 


0.42±0.03 


0.18±0.01 


0.18±0.01 


0.18±0.01 


[Sll]A4076 










0.091±0.006 


0.17±0.01 


0.18±0.01 


0.050±0.004 


0.054±0.004 


H(5A4102 


0.31±0.02 


0.24±0.02 


0.34±0.02 


0.24±0.02 


0.23±0.02 


0.35±0.03 


0.19±0.01 


0.18±0.01 


0.18±0.01 


H7A434O 


0.47±0.03 


0.44±0.03 


0.56±0.04 


0.41±0.03 


0.45±0.03 


0.64±0.05 


0.40±0.03 


0.43±0.03 


0.41±0.03 


[Oiii]A4363 










0.099±0.007 


0.14±0.01 


0.090±0.006 


0.092±0.006 


O.lliO.Ol 


[Oiii]A4959 


0.88±0.06 


1.11±0.08 


0.92±0.07 


1.38±0.10 


1.30±0.09 


1.54±0.11 


1.30±0.09 


1.23±0.09 


1.48±0.10 


[Oiii]A5007 


2.54±0.18 


3.25±0.23 


2.74±0.19 


4.64±0.33 


3.73±0.26 


4.55±0.32 


3.61±0.26 


3.53±0.25 


4.22±0.30 


F(H/3) 


1.17±0.06 


2.61±0.13 


1.72±0.09 


2.18±0.11 


3.74±0.19 


1.32±0.07 


2.43±0.12 


4.28±0.21 


6.84±0.34 


EW(H/3) 


52±3 


88±4 


156±8 


98±5 


50±3 


32±2 


72±4 


72±4 


89±4 


EW([Oll] A3727) 


109±5 


117±6 


151±8 


112±6 


88±4 


80±4 


88±4 


77±4 


77±4 


EW([Oiii] A4959) 


43±2 


78±4 


417±21 


136±7 


60±3 


67±3 


71±4 


96±5 


148±8 


EWQOiii] A5007) 


97±5 


194±10 


437±20 


400±20 


165±8 


189±10 


149±7 


219±11 


322±16 



H/3 flux in 10 erg cm ^ sec ^; EW in A. 



Table 3. Observed "Red" Emission Line Intensities Relative to Ha. 



Galaxy 


F563-1 


F563-1 


F563-1 


F571-5 


U1230 


U1230 


U1230 


U1230 


U1230 


Hii 


1 


2 


3 


1 


1 


2 


3 


5 


6 


He IA5876 




0.043±0.003 








0.053±0.004 


0.055±0.004 






[Ol]A6300 




















[Nll]A6548 


0.024±0.002 


O.OlOiO.OOl 


0.033±0.002 


0.045±0.003 




0.022±0.002 


0.022±0.002 






[Nll]A6584 


0.082±0.006 


0.055±0.004 


0.069±0.005 


0.086±0.006 


0.052±0.004 


0.087±0.006 


O.lliO.Ol 


0.092±0.007 


0.12±0.01 


[S Il]A6717 


0.17±0.01 


0.081±0.006 


0.19±0.01 


0.053±0.004 


O.lSiO.Ol 


0.19±0.01 


O.lSiO.Ol 


0.19±0.01 


0.13±0.01 


[S ii]A6731 


0.12±0.01 


0.072±0.005 


0.13±0.01 


0.025±0.002 


0.14±0.01 


0.14±0.01 


O.lliO.Ol 


0.13±0.01 


0.17±0.01 


He IA7065 




0.019±0.001 
















F(Hq) 


2.67±0.13 


4.45±0.22 


0.95±0.05 


3.27±0.16 


0.76±0.04 


3.10±0.16 


4.92±0.25 


0.98±0.05 


0.53±0.03 



Ha flux is in 10 erg cm ^ s 



Table 3 - continued 



Galaxy 


U5005 


U5005 


U9024 


U9024 


U12695(l) 


1112695(1) 


U12695(2) 


U12695(2) 


U12695(2) 


Hii 


2 


3 


1 


3 


2 


3 


2 


3 


4 


He IA5876 


0.032±0.002 


0.034i0.002 




0.043i0.003 


0.041i0.003 


0.057i0.004 


0.040i0.003 


0.034i0.002 


0.036i0.003 


[Oi]A6300 


0.031±0.002 


0.025i0.002 






0.017i0.001 




0.026i0.002 


0.033i0.002 


O.OlOiO.OOl 


[Nll]A6548 




















[Nll]A6584 


0.066±0.005 


0.056i0.004 


0.039i0.003 


0.061i0.004 


0.034i0.002 


0.032i0.002 


0.051i0.004 


0.030i0.002 


0.015i0.001 


[S ll]A6717 


O.lSiO.Ol 


O.lliO.Ol 


0.14i0.01 


0.082i0.006 


0.083i0.006 


0.094i0.007 


0.12i0.01 


0.080i0.006 


0.041i0.003 


[S Il]A6731 


O.lliO.Ol 


0.076i0.005 


0.092i0.007 


0.063i0.004 


0.034i0.002 


0.045i0.003 


0.038i0.003 


0.033i0.002 


0.017i0.001 


He IA7065 




0.013i0.001 






0.012i0.001 




0.042i0.003 


0.017i0.001 




F(Ha) 


3.26i0.16 


7.53i0.38 


2.77i0.14 


5.20i0.26 


ll.51i0.58 


3.17i0.16 


7.72i0.39 


ll.08i0.55 


13.15i0.66 



Ha flux is in 10 erg cm s ^ 
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Figure 8. Diagnostic diagram for the H II regions of tlie cuiTent sample. 
H II regions sliould fall along the line from Baldwin, Phillips & Terlevich 
(1981). 



Figure 9. Comparison of oxygen line ratios for the current sample of H II 
regions to the observations of McGaugh (1994). To separate the [Oil] and 
[O III] data, a shift of 0.05 dex has been applied to the x-axis for the [O III] 
data. 



H II regions have been observed on both occasions. The two 
outUers are UGC 1230 #6 and UGC 12695(1) #3. UGC 1230 #6 
is a faint H II region and the uncertainties in the original McGaugh 
(1994) data are large. The current measurement of the H/3 flux is 
much less than the measurement of McGaugh (1994). The slit may 
have missed the majority of this H II region during the more recent 
observation. All of the H II regions for both slit positions and both 
observations of UGC 12695 are large and bright. Because the H II 
regions are large, it is possible that they may resolve into smaller 
ones. Either the slit fell across a different H II region entirely, or a 
different part of the same H II region during the two observations. 
For these reasons, we exclude Hll regions UGC1230 #6 and UGC 
12695(1) #3 from further analysis. Together, Figs. 8 & 9 demon- 
strate for the rest of the H II regions that even though we cannot 
intercompare lines from red and blue spectra, we can construct 
valid line ratios from within each. 



4 OXYGEN ABUNDANCE 

4.1 [O III] A4363 Direct Abundance Measure 

Standard methods of determining oxygen abundances involve 
measuring a temperature sensitive line such as [O III] A4363 (Os- 
terbrock 1989). Ionic abundances of O"''^ and can then be de- 
rived using the electron temperature. We have only 4 H II regions 
which have both an [O III] A4363 measurement and a reddening 
correction available from McGaugh (1994) (see Section 4.2). We 
determine the temperature using the low-density limit (as indicated 
by the density-sensitive [S II] A6717/6731 ratio) of the relation 
between the [Olll] (A4959 -I- A5007)/A4363 ratio and the tem- 
perature (Osterbrock 1989). We then follow the method of Pagel 
et al. (1992) to determine the total oxygen abundance using both 
0++/H+ and 0+/H+. The abundances based on the [O III] A4363 
line are listed in Table 4. Despite being applicable to only 4 re- 
gions, the abundances determined in this fashion can provide an 
interesting check on abundances derived from indirect techniques. 



Table 4. Direct Oxygen Abundance Measurements 



Galaxy & region 


log (O/H) 


UGC 12695(1) 2 


-4.27±0.07 


UGC 12695(2) 2 


-4.30±0.06 


UGC 12695(2) 3 


-4.35±0.06 


UGC 12695(2) 4 


-4.32±0.04 



4.2 i?23 strong Line Method 



[O III] A4363 is notoriously faint, and for much of the present data 
set is undetected. When this is the case, one must rely on alternate 
methods of determining abundances. A widely accepted approach 
is the "R23" strong line method. As suggested by Pagel e t alj 
J1979'). the oxygen abundance (O/H) can be determined from the 
sum of the [O II] A3727, [O III] A4959 and [O III] A5007 intensities 
relative to H/3. 

The oxygen line strength forms two distinct branches. At 
higher abundances, cooling from the infrared lines is important 
and the sum of the optical [O II] and [O III] lines is low. As the 
abundance decreases, cooling comes increasingly from the [O II] 
A3727 and [O III] AA4959,5007 lines, and the sum of these lines 
increases. This continues until oxygen becomes rare, and at this 
point the sum of the [O II] and [O III] lines begins to once again 
decrease. The degeneracy between these branches is broken by the 
ratio of [N II] A6584 to [O II] A3727, although other Une pairs may 
be used (eg. Alloin et al. 1979; van Zee et al. 1998). R23 is: 



R2i 



[O II] A3727 + [O III] A4959 -f [O iii]A5007 

Wp ■ 

Also frequently used in conjunction with R23 is O32: 
[Oiii]A4959 + [Oiii]A5007 



O32 



[Oii]A3727 



(1) 



(2) 



an indicator of the ionization parameter. 

There have been several proposed calibrations of the R23 
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method; we use the calibration of'McGaueh' ('l99l') as represented 
by the analytic functions described in Appendix A. These return 
the oxygen abundance and ionization parameters once the branch 
of the R2S calibration is specified by the ratio of [N Il]/[0 II]. The 
pure model fit has been used in the following analysis. 

Due to the limitations imposed by the inability to link to- 
gether the blue and red spectra, we encountered two problems 
when applying the R23 strong line method. First, we cannot sim- 
ply use the ratio of [N II] A6584 to [O II] A3727 to break the branch 
degeneracy-the lines are on separate arms of the spectrum. Instead 
we begin by using the [Nil] A6584 to Ha ratio (van Zee et al. 
1998). Crudely speaking, a ratio greater than 0.1 indicates the up- 
per branch and less than 0.1, the lower branch. Then, ignoring our 
calibration difficulties, we determine the branch with the [N II] to 
[O II] ratio and compare the results. In all cases except UGC 1230 
#5, the branch determination is the same regardless of the method 
used. The results using the [Nil] to Ha ratio are also the same 
when compared to McGaugh (1994). The exception to this is UGC 
1230 #3 which is designated as a lower branch object by McGaugh 
(1994) but is determined here to be on the upper branch. For fur- 
ther analysis, we adopt the lower branch designation for both UGC 
1230 #3 and #5. The second problem encountered is that the R23 
strong line method requires that the line fluxes be corrected for 
reddening. We are unable to use the Ha/H/3 ratio to determine red- 
dening corrections because the two lines fall on separate arms of 
the spectrum. We attempted to use H/3/H7 to obtain the redden- 
ing, but the uncertainties were very large. We therefore adopted 
the reddening corrections from McGaugh (1994). This enabled the 
R23 strong line method to be applied to 9 H II regions. The results 
are listed in Table 5. 



4.3 Equivalent Width Method 

There are 7 Hll regions for which a reddening correction from 
McGaugh (1994) is unavailable. For these regions (as well as 
the others of the current sample), we use an equivalent width 
version of the R23 method to obtain the abundance. Kobulnicky 
& Phillips (2003) found that in the absence of a reddening cor- 
rection, equivalent width (EW) ratios are an adequate substi- 
tute for reddening-corrected flux ratios of the [Oil], [Olll], and 
H/3 lines in the R23 method. To correct for stellar absorption, 
a 2 A correction was made to the EW of H/3 for each Hll 
region jM cCall. Rvbski, & Shields 1985; Oe v & Kennicullll993l: 
lMcGaugh.l994;.Kobulnickv & PhiniDS..2003i) . Following Kobul- 
nicky & Phillips (2003), the R23 strong line method was then im- 
plemented using the measured equivalent widths of the lines rather 
than the fluxes. Equations 3 and 4 are the equivalent width versions 
of R23 and O32 : 

rpTJ/ - -^'^[OlIlA3727 + ^^[O III1A4959 + -^^[O III]A5007 
J^Wr23 = ;:;777 (3) 



EWo, 



EWhp 

-^^[OlII]A4959 + -^^[OlII]A5007 
^"^{O II]A3727 



(4) 



The branch was determined using the ratio of the observed, un- 
corrected [N II] A6584 to Ha flux. The equivalent widths of H/? 
and the [O II] and [O III] lines are listed in Table 2 and the result- 
ing oxygen abundances are in Table 5. We also apply the equiva- 
lent width method to a number of other H II regions of McGaugh 
(1994); the equivalent widths and corresponding oxygen abun- 
dances of those regions are listed in Appendix B. 



Table 5. Oxygen 


Abundances from Indirect Methods 




Galaxy & region 


log(0/H) 


log(0/H) 


loH(0/H) 




(EW) 


(R23) 


(McGaugh 1994) 


F563-1 1 


-3.87±0.05 






F563-1 2 


-3.70±0.06 






F563-1 3 


-3.53±0.06 






F571-5 1 


-4.11±0.05 






UGC 1230 I 


-4.49±0.04 


-4.21±0.10 


-4.17±0.30 


UGC 1230 2 


-4.07±0.05 


-4.11±0.07 


-4.04±0.16 


UGC 1230 3 


-4.20±0.04 


-4.17±0.03 


-4.16±0.09 


UGC 1230 5 


-4.20±0.05 


-4.12±0.I5 


-3.66±0.47 


UGC 5005 2 


-4.22±0.04 






UGC 5005 3 


-4.35±0.04 






UGC 9024 1 


—4 90+0 05 




-3.95±0.21* 


UGC 9024 3 


-4.17±0.05 


-3.79±0.23 


-3.70±0.36 


UOL- 12oVj(1 } L 


-4.13±0.05 


-3.93±0.04 


A no_l_A no 


UGC 12695(2) 2 


-4.40±0.04 


-4.04±0.09 


-4.01±0.12 


UGC 12695(2) 3 


-4.30±0.05 


-4.12±0.04 


-4.11±0.08 


UGC 12695(2) 4 


-4.26±0.05 


-4.13±0.06 


-4.14±0.08 



Column I abundances are determined using the equivalent width method. 
Column 2 abundances use the i?23 method & reddening connections 
adopted from McGaugh (1994). Column 3 abundances are from McGaugh 
(1994). 

Errors are due to propagation of measurement uncertainties and do not in- 
clude calibration uncertainties of at least 0.1 dex. * Ambiguous branch 
determination; lower branch abundance is listed here. 

5 ABUNDANCE RESULTS AND DISCUSSION 
5.1 H II Region Abundances 

We now have at least 1 and up to 4 measurements of the oxy- 
gen abundance for each H II region. Only the quality of the data 
should affect the abundances using the ii23 method, as both the 
method and calibration are identical. The EW method and the 
[O III] A4363 direct method, however, are distinct and those results 
should be compared to the R23 method. The abundances derived 
from the [O III] A4363 measurements agree with the EW method 
abundances within the uncertainties. As was found by McGaugh 
(1994), the direct oxygen abundances are somewhat lower than 
those of the R23 method. In Fig. 10 we plot the 7?23 -derived abun- 
dances of the current data and the McGaugh (1994) data against the 
EW-derived abundances. In almost all cases the 7?23 abundances 
are higher than the EW abundances. The abundances of McGaugh 
(1994) are higher than the abundances of the current data, but agree 
to within the errors. 



5.2 Mean LSB Galaxy Abundances 

Ideally, we would like to have enough information to determine 
abundance gradients and use them to establish a characteristic 
metallicity for each of our galaxies. This, however, is a difficult 
challenge in LSB galaxies. Unlike high surface brightness galax- 
ies where bright H II regions are plentiful, most LSB galaxies lack 
enough bright H II regions to make an abundance gradient deter- 
mination (Figs. I — 6). We therefore choose to take the mean of 
the abundances of the individual H II regions as representative of 
the entire galaxy. One should bear in mind that this mean value 
may or may not be a good indicator of the true global metallicity. 
Confidence in the use of the mean abundance as the overall galaxy 
abundance seems to be appropriate, though, as de Blok & van der 
Hulst (1998) find no abundance gradient, but instead a constant 
abundance with radius, for 3 LSB galaxies which had enough H II 
regions distributed over the length of the galaxy to make a gradi- 
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Figure 10. Comparison of oxygen abundances derived using the R23 
metliod and tlie EW method. Lines connect the regions which are part of 
both the current sample and the McGaugh (1994) sample. Crossed lines in 
the lower right represent typical error bars. 

ent study feasible. It should also be noted that abundances correlate 
with local surface density as well as with radius (McCall 1982). By 
definition, LSB galaxies have lower surface densities and little dy- 
namic range in the surface brightness over which H II regions are 
observed. One would thus expect little change in metallicity across 
an LSB galaxy. With this in mind, we take the average of the H II 
region abundances in each galaxy to represent the global abun- 
dance of that galaxy. Table 6 lists both the average EW abundance 
and the average R23 abudance. The data for our current sample 
of LSB galaxies has been supplemented with data from McGaugh 
(1994) and de Blok & van der Hulst (1998). The McGaugh (1994) 
data are for 3 H II regions from our galaxy sample that were not ob- 
served as part of the current sample. The de Blok & van der Hulst 
(1998) data are R23 method oxygen abundances for different Hll 
regions of many of the same galaxies as investigated here. 



6 LUMINOSITY-METALLICITY RELATION 

The luminosity-metallicity (L-Z) relation, and mass-metallicity re- 
lation, have been important tools for studying galactic chemical 
evolution. Relations have been found for, among others, dwarf Ir- 
regulars (eg. Lee et al. 2003), late-type galaxies (eg. Pilyugin & 
Ferrini 2000) and ellipticals (eg. Brodie & Huchra 1991). L-Z re- 
lations have also been found using very large survey samples (eg. 
KISS Melbourne & Salzer 2002; 2dFGRS Lamareille et al. 2004; 
SDSS Tremonti et al. 2004). The correlation between luminosity 
and metallicity exists over a wide range of luminosities and metal- 
licities, although there are indications that perhaps dwarf galaxies 
have a shallower L-Z relation and that one L-Z relation is there- 
fore insufficient to describe all galaxy types. It is interesting to see 
where LSB galaxies fit in. 

Fig. 1 1 is the L-Z plot for the LSB galaxies. Also shown are 
a variety of recent fits to the L-Z relation. There is some disparity 
here which may stem in part from the use of different methods 
and calibrations for estimating abundances. It appears that the LSB 



Figure 11. Luminosity-Metallicity plot for LSB galaxies. The thin solid 
line is from the KISS sample of Melbourne & Salzer (2002). The heavy 
solid line is from the 2dFGRS sample of Lamareille et al. (2004). The 
short-dash line is the Tremonti et al. (2004) SDSS sample, and the long- 
dash line is for the dirr sample of Lee et al. (2003). The majority of data 
from the Tremonti et al. (2004) fit have Mg brighter than —18. The fit is 
extrapolated to fainter magnitudes in their fig. 4, and further extrapolated 
here. 



data prefer a steeper relation similar to that describing the 2dFGRS 
sample. 

Because it is more difficult to estimate mass, mass-metallicity 
relations are less often computed than L-Z relations. We have mass 
estimates for our sample, and present the mass-metallicity plot in 
Fig. 12. Also plotted in Fig. 12 is the dIrr sample of Lee et al. 
(2003) and the results of the mass-metallicity relation determined 
for SDSS galaxies of Tremonti et al. (2004). Tremonti et al. (2004) 
find the correlation between mass and metallicity to be roughly 
linear between 10* '^ Mq and 10^" '^ Mq before it begins to flatten 
out. The LSB galaxies mainly fall within this linear region, but 
have considerably lower abundances in the mean. 

The calibration of both mass and metallicity is slightly differ- 
ent here than in Tremonti et al. (2004). The R23 calibration used 
here would give lower (O/H), while the IMF used here is 0.15 dex 
more massive (Bell et al. 2003). This would shift the SDSS line 
down and to the right in Fig. 12. While the shift in mass is trivial 
to achieve, that in abundance is more complex because the shape of 
the calibrations are somewhat different. For now we ignore these 
differences as the net difference is not large compared to the range 
of the plot. One should bear in mind the caveat that it is impossi- 
ble to exclude the possibility that the apparent differences might 
be reconciled by a change in calibration (see appendix A). 

Accepting the calibration of the mass-metalicity relation as 
published by Tremonti et al. (2004), the mean abundance of both 
LSB galaxies and the dIrr sample of Lee et al. (2003) is lower 
than that of the SDSS galaxies. This is qualitatively consistent with 
higher gas fraction galaxies being less evolved at a given current 
stellar mass. Certainly it is possible that LSB galaxies simply re- 
side on the low metallicity side of the overall distribution, a result 
which is consistent with previous work (McGaugh 1994; de Blok 
& van der Hulst 1998). However, the mass-metallicity relation is 
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Table 6. Mean Oxygen Abundances and Gas Mass Fractions 



Galaxy <log(0/H)> <log(0/H)> Ref. Mb (B-V) MriILb fg Method 
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References for mean R23-derived abundances— 1: CuiTent data; 2: McGaugh (1994); 3: de Blok & van der 
Hulst (1998) 

Gas Mass Fraction Methods— 4; de Blok & McGaugh (1998); 5: Uses Eq. 7 of the text and relevant data 
from McGaugh & de Blok (1997); 6: McGaugh & de Blok (1997) 
* Mhi/Lr from van der Hulst et al. (1993) 
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Figure 12. Mass-Metallicity plot for LSB galaxies. LSB galaxies are rep- 
resented by circles. The stars are the dirr sample from Lee et al. (2003) 
and the contours are the results of Tremonti et al. (2004) for 53,000 SDSS 
star-forming galaxies. The contours enclose 68% and 95% of their data. 



not well defined below log Af, = 8.5, so it is also possible that the 
relation turns down sharply at this point. 



7 OXYGEN ABUNDANCE VS. GAS MASS FRACTION 

Additional clues to galactic chemical evolution can be gleaned 
from the gas mass fraction. The gas mass fraction (fg) is a mea- 



sure of the total mass in gas (Alg — molecular, atomic, ionized and 
metals) to the combined total mass in gas and stars (M*): 



Ma 



M„ + M* 



Rephrased in terms of observables, the gas mass fractions is 



(5) 



(6) 



where T« is the stellar mass-to-light ratio (M«=T«L), L is the 
luminosity, AIhi is the mass in neutral hydrogen and rj serves as 
the equivalent to T« {Mg=riMHi)- 

Together, metallicity and the gas mass fraction can be used as 
a diagnostic for testing different models of the chemical evolution 
of galaxies. In the closed-box model, there is neither inflow nor 
outflow of mass as the initially metal-free gas evolves with time. 
Mathematically, this is represented as 



(7) 



where Z is the metallicity and p is the yield (eg. Searle & Sargen t 
1972; Tinsley 1980; Edmunds 1990). Following lLee et alj iioosi) . 
this can be rewritten in a more user-friendly fashion: 



12 + log(0/H) = 12 + log(0.196yo) + log log(l/fg 



(8) 



where yo is the oxygen yield. Written this way, the closed-box 
model predicts a slope of unity; deviations from this slope may be 
an indication for the inflow or outflow of gas. 

Our best estimates of the gas fractions for the entire disk are 
listed in Table 6. We give preference to stellar mass estimates from 
de Blok & McGaugh (1998). For galaxies not in that sample, we 
compute the gas fractions with Eq. 6 using 7^=1.4 and compute T* 
from {B — V^) using 



log(A/./Ls) = 1.737(5 - V) - 0.787 



(9) 
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Figure 13. Plot of gas mass fraction versus metallicity. Circles are the cur- 
rent data s et of L SB galaxies. Stars are data for field dwarf Irregulars from 
iLee et all i2003 V The solid line indicates the behaviour of the closed box 
model with a yield of 0.002. 

iBell et alfcOO Sh (their 'diet' Salpeter IMF). We choose this IMF 
because it is fairly close to optimal with respect to dynamical data 
(McGaugh 2004). Finally, if neither of the two previous methods 
can be applied, we assume the gas fraction in table 1 of McGaugh 
& de Blok (1997). Plotted in Fig. 13 is oxygen abun dance ver- 
sus ga s fraction for our data and data for field dirr from lLee et alj 
i2003h . There is no clear evidence of deviation from the closed- 
b ox model evident in t his plot. Our data are consistent with those 
of lLee et al1i2003h and lPagell fl986 ) albeit with more scatter. 

LSB galaxies are consistent with the yield of p=0.002 deter- 
mined for other galaxy types (Lee et al. 2003; Pagel 1986). This 
is an interesting result for the LSB galaxies as it is often argued 
that they hav e shallow poten tial wells, making it easier for mass 
loss to occur. iGarnetll i2002ft finds that outflows may be impor- 
tant in galaxies with Vrot<50 km s^^. All galaxies discussed here 
have Vrot>50 km s~^. The apparent uniformity of the yield is also 
suggestive of a universal IMF, as favored by dynamical constraints 
(McGaugh 2004). However, the abundance data alone leave sub- 
stantial room to consider other models, such as the bottom-heavy 
IMF for LSB galaxies suggested by Lee et al. (2004). 



8 CONCLUSIONS 

We have studied the optical emission line spectra of H II regions 
in six LSB galaxies. We have supplemented this sample with the 
LSB galaxy data of McGaugh (1994) and de Blok & van der Hulst 
(1998). We have measured oxygen abundances using both the EW 
method and R23 strong line method, and in a few cases directly us- 
ing the [O III] A4363 temperature sensitive line. The oxygen abun- 
dances of the LSB galaxies show them to be metal poor, while 
the gas mass fractions are high. High gas mass fractions and low 
metallicities are consistent with LSB galaxies being at an early 
stage of their evolution. There is no need to suppose that they 
are fundamentally different from other types of galaxies. In gross 
terms, a closed-box model with the same effective yield as for other 
galaxies appears to be adequate. There is no evidence in the gas 



fractions that requires extensive gas loss stemming from blow-out 
in these systems. Finally, we have found that LSB galaxies are 
broadly consistent with the luminosity-metallicity relation of other 
galaxy types. 
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APPENDIX A: ANALYTIC EXPRESSIONS FOR STRONG 
LINE OXYGEN ABUNDANCES 

The i?23 method <PageletalJll979h is a valuable tool for obtain- 
ing oxygen abundance estimates. The strong lines upon which it 
is based are readily accessible observationally. While the temper- 
ature sensitive [O III] A4363 line can usually only be observed in 
bright, nearby, low metallicity H II regions, the strong lines are of- 
ten measurable to cosmologically significant distances. 

To facilitate the use of the R2S method for large data sets, 
we provide here sets of equations relating the observables _R23 and 
O32 to the oxygen abundance and ionization parameter. To sim- 
plify the appearance of these relations, we employ the following 



notation: 

Z = log(0/H) (Al) 

U = logU (A2) 

X = log(i?23) (A3) 

y = log(032). (A4) 



Z is used here as a shorthand for the logarithmic oxygen abun- 
dance, and should not be confused with the mass fraction of metals. 
U is the volume averaged ionization parameter, being basically the 
ratio of the number of ionizing photons to the number of particles. 

Below we provide functions Z[x,y) and U{x,y) that give 
the nebular parameters. Two similar but distinct calibrations are 
provided. The first is a straight representation of the model grid 
of lMcGaughl Il99ll) . The second attempts a modest correction to 
this grid to account for the observed decrease in stellar effective 
temperatures at high metallicity by fitting to the metal rich H II re- 
gion Searle 5 ^Shields & Searle* 1978; Kinkel & Rosa 1994). These 
relations have been employe d to e stimate oxygen abun danc es in 
high redshift galaxies b y iLillv. Carollo. & StocktonI i2003l) and 
iKobulnickv e'tai]<2003h . 

Al Model Calibration 

A functional repr esentation of the Mu = 60 A/q model grid of 
iMcGaugi] il99lh is given here. The i?23-0/H relation is double 
valued, so we split it into separate expressions for the upper (high- 
Z) and lower (low-Z) branches. On the lower branch. 



Zi{x, y) = -4.93 + 4.25a; - 3.35 sin{x) 

-0.26i/-0.12sin(i/) (A5) 

Ui{x, y) = -2.95 + 0.17x2 + 1.022/, (A6) 

while on the upper branch, 

Z„{x,y)= -2.65 -0.91X + 0.12?/ sin(x) (A7) 

!7„(x, y) = -2.39 - 0.35x + 1.29y - 0.15a:y. (A8) 



These expressions reproduce the model grid to within 0.1 dex 
when the arguments of the trigonometric functions are treated as 
being in radians. 



Distinguishing which branch is appropriate can be accom- 
plished with the [N II] A6583 line iMcGaugh 199lLll994 . Strong 
(weak) nitrogen emission indicates that the upper (lower) branch 
is appropriate. The dividing line is at log([N n]/[0 II]) « — 1. 
This can often be judged by eye in raw spectra: is [N II] prominent 
compared to Ha, or is it very weak? However, there is no precise 
quantitative dividing line, as the [N II]/Ha ratio depends on a va- 
riety of factors. 

The [N II]/[0 II] ratio is itself a very good abundance indica- 
tor. In two respects it is superior to i?23: (1) it is monotonic rather 
than double- valued, and (2), being formed from lines of the same 
ionization state, is insensitive to the ionization parameter U. Un- 
fortunately it is rather difficult to calibrate because of the behavior 
of N/O at low O/H, and in practice is very sensitive to the red- 
dening correction because of the large separation of the lines in 
wavelength. 

A2 Semi-Empirical Calibration 

A large variety of _R23 calibrations have been suggested since 
the met hod was introduced by Pagel (e.g, see the compila- 
tions of iMcGauehl Il99ll : IZaritskv. Kennicutt. & Huchrj Il994l: 
iKewlev ^^ODita200 ?i~Some of these calibrations are empirical, 
some are based on models, and others are merely hybrids of previ- 
ous calibrations. The plain fact is that it remains difficult to pin an 
absolute calibration on i?23. 

The problem of absolute calibration is particularly acute at 
high (> solar) metallicity where there is a dearth of independent 
data which test the method. On the one hand, it is clear that R23 
is very sensitive to oxygen abundance, and hence a good relative 
abundance indicator. At and above solar metallicity, slight changes 
in O/H cause large ones in i?23, so the discrimination is very pow- 
erful. However, the power to discriminate between slightly differ- 
ent oxygen abundances does not extend to knowledge of the abso- 
lute abundance. 

We provide here a semi-empirical calibration which attempts 
to fit the one well kno wn high metallicity p oint Searle 5. The orig- 
inal determination ( Shields c& Searle 1978) of the abundance of 
Searle 5 is Z = -2.89. Kinkel & Rosa ( 1994) estimate a lower 
oxygen abundance of Z = —3.12. T he latter deter mination is con- 
sistent with the detailed modeling o f lEvanslJ 19861) . 

Nearly all R23 calibrations in the literature overestimate 
the lower abundance measurement, including the model grid of 
iMcGaughl tl99D as given in the previous section. The model grid 
of lMcGaughl i 199 ill does give this lower abundance if one takes 
« 38, OmK as inferred bv lShields & Searld <1978h and sub- 
sequent workers. The chief problem seems to be in the ionizing 
spectra of metal rich stars which are needed as input to the pho- 
toionization models. This is a highly nonlinear problem beyond 
the scope of this paper. Nevertheless, we can interpolate through 
the existing model grid to follow the observed drop in with Z. 
In effect, we have added an additional term to the T-t,[M, Z) rela- 
tion discussed bv lMcGaughl (19^|) to empirically compensate for 
the excess hardness of the non-LTE stellar atmospheres employed 
in those models. 

The result is a semi-empirical calibration which splits the dif- 
ference between the Shields & Searle (1978) and Kinkel & Rosd 
(1994) measurements, yielding Z = —3.03 for Searle 5. This is 
within 0.1 dex of the iKinkel & Rosal il994i[) measurement, which 
is about what we expect for the accuracy of this approach. Polyno- 
mial fits to the semi-empirical calibration are 

Ze{x, y) = -4.944 + 0.767x + 0.602a;2 

-y(0.29 + 0.332a; - 0.331a;2) (A9) 
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and 

Zu{x, y) = -2.939 - 0.2a: - 0.237x^ - 0.305x^ 

-0.0283X'' - y(0.0047 - 0.0221a; (AlO) 
-0.102a;2 - 0.0817a;^ - 0.00717a;'*). 

The terms involving xy and x'^y in the expression for Zu{x,y) 
have little impact on the fit but are included for completeness. 

Note that U and its indicator O32 play an important second 
parameter role on the lower branch which persists to rather high 
abundances (Z ~ Zq, x ~ 0.5). This effect is not suppressed. 
In fact, it is particularly important around Z ~ 1/2Zq, x > 0.7. 
Many data lie in this region where other calibrations rely on only a 
single parameter. However, at still higher abundances, the sensitiv- 
ity of R23 to the oxygen abundance becomes so great that it effec- 
tively becomes a one parameter sequence. If there is any residual 
second parameter dependence on the upper branch, it is more likely 
due to T^, than U. 

The net result of this modest recalibration of the model grid 
is a slower variation of the oxygen abundance with R23 rela- 
tive to most other calibrations. For example, weak R23 lines with 
a; = give Z — —2.9 3, 0.2 dex less than the calibrations of 
lEdmunds & PagelHl984 . McCa ll et al. (1985) and Zari tsky et al. 
(1994), but comparable to that of iDopita & Evans! jl98d) . 

The expressions given here provide a uniform, self-consistent, 
accurate relative calibration of the R2S abundance sequence. It is 
difficult to estimate the error in the absolute calibration: we do 
not know what the intrinsic i?23-0/H relation actually is. Our cal- 
ibration does at least provide a consistent, easy to use method for 
estimating abundances over a very large range of parameters, ex- 
plicitly incorporating the important effects of ionization. 

An additional note: Long after this appendix was written, 
and very late in the refereeing process, a paper appeared containing 
more direct estimates of the oxygen abundance in high abundance 
Hll regions in M51 (Bresolin, Garnett, & Kennicutt 2004). The 
semi-empirical calibration above, while lower than most other cal- 
ibrations and consistent with the one previously well known case, 
overestimates the abundances of Bresolin et al. (2004) by ~ 0.2 
dex. The R23 calibrations discussed here are still valid indicators 
of the relative oxygen abundance, but the results of Bresolin et al. 
(2004) highlight the difficulty of obtaining an absolute calibration. 
Since the strong line method also tends to overstate lower branch 
abundances by a comparable amount, a simple shift might be 
crudely acceptable. This would amplify the disparity between the 
LSB galaxy abundances determined here and the mass-metallicity 
relation of Tremonti et al. (2004). If the results of Bresolin et al. 
(2004) are correct, the calibration used by Tremonti et al. (2004) 
may overstate abundances by as much as 0.5 dex. This might go 
some way to explaining why their mass-metallicity relation implies 
super-solar abundances for Milky Way mass galaxies. Further in- 
vestigation of the calibration of the strong line method, especially 
along the upper branch, is clearly warranted. 



APPENDIX B: EQUIVALENT WIDTH MEASUREMENTS 
OF R23 LINES OF MCGAUGH 1994 

In Table BI we list the equivalent widths of the H/3, [Oil] and 
[Qui] lines of the Hll regions of McGaugh (1994) which supple- 
ment the current sample of LSB galaxies. The equivalent widths 
are used to determine the oxygen abundance, also listed in Ta- 
ble Bl, using the EW method described in Section 4.3. For com- 
parison, the McGaugh (1994) oxygen abundances using the R23 
method are also listed. 



Chemical Evolution ofLSB Galaxies 

Table Bl. Equivalent Widths and Oxygen Abundances of McGaugh (1994) 



Galaxy & region H/3 [Oil] [OlII] [OlII] log(0/H) log(0/H) 
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Equivalent widths in A. 

'Ambiguous branch determination; lower branch abundance is listed here. 
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